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A  unit  vector 
Any  continuous  function 
Enthalpy  per  unit  mass  of  fluid 
Total  enthalpy  per  unit  mass  of  fluid 
Relative  total  specific  rothalpy 

Total  enthalpy  of  a  fluid;  in  linear  systems,  the  total  enthalpy  of  a  fluid 
per  unit  length  of  blade 

Lift  per  unit  length  of  blade 

Mass  flow  rate 

A  constant 

Local  static  pressure 

Total  pressure 

A  generalized  curvilinear  coordinate 
A  specific  heat  rate 
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ABSTRACT 


A  differential  turbomachinery  equation  describing  the  energy  transfer  between  o 
fluid  and  any  body  moving  in  that  fluid  was  derived.  The  derivation  is  based  upon  the 
Coriolis  form  of  the  Navier-Stokes  Equations.  A  differential  equation  for  the  total  relative 
rothalpy  is  also  obtained.  Both  equalities  contain  a  rigorous  viscous  correction.  Both  may 
be  evaluated  in  the  absolute  and  moving  frame. 

On  integration  of  the  differential  equations,  a  form  of  the  Euler  Turbomachinery 
Equation  with  viscous  correction  is  derived.  The  resultant  form  contains  two  distinct  work 
rate  terms  for  the  axial  and  radial  components  of  the  flow.  The  fact  that  integration  yields 
a  result  which  approximates  the  classic  Euler  Turbomachinery  Equation  constitutes  con¬ 
firmation  of  the  derivation. 

An  application  of  the  eauation  to  an  ideal  infinite  linear  cylinder  with  bound  vortic- 
ity  was  developed,  yielding  the  expected  known  result. 

ADMINISTRATIVE  INFORMATION 

This  study  is  a  further  development  of  earlier  work  (1987)  performed  in  partial  ful¬ 
fillment  of  thesis  requirements  of  the  Aerospace  Department  of  the  University  of 
Maryland  at  College  Park.  Professor  Everett  Jones  was  the  graduate  advisor  and  director. 
Publication  of  this  manuscript  was  supported  by  the  David  Taylor  Research  Center  under 
the  Independent  Research  and  Development  Program  (IR&IED)  of  Dr.  Bruce  Douglas. 

The  work  has  significant  implications  to  turbomachinery  studies  of  efficiency  and 
noise  in  the  Power  Systems  Division,  Code  272,  of  the  Propulsion  and  Auxiliary  Systems 
Department. 
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INTRODUCTION 

A  monument  of  tuibomachinery  technology  is  the  Turbomachinery  Equation  (ulti¬ 
mately  due  to  Euler’s  analysis  of  torque  in  fluids)  which  is  based  upon  thermodynamic 
definitions  of  work  and  Newton’s  Laws.  Since  the  Navier-Stokes  Equations  and  Crocco’s 
Equation!  1,2]  in  a  rotating  (moving)  frame  are  also  based  upon  thermodynamics  and 
Newton’s  Laws,  they  must  in  principle  contain  the  Turbomachinery  Equation  in  differen¬ 
tial  form  and,  on  integration,  in  integral  form. 

The  Coriolis’  transformation  relates  moving  rotor  frames  and  absolute  or  laboratory 
frames  so  that  integration  of  the  energy  rate  may  be  conveniently  performed  in  a  time-in¬ 
dependent  frame  with  a  time-independent  set  of  coordinates. 

The  transformation  leads  to  simplified  expressions  for  the  substantial  derivative  of 
the  total  enthalpy  uncoupled  from  the  substantial  derivative  of  the  entropic  energy.  The 
uncoupling  generates  a  differential  Turbomachinery  Equation  which,  on  integration, 
yields  a  novel  form  of  the  Turbomachinery  Equation  corrected  for  viscous  non-ideal  flow. 
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The  classical  Turbomachinery  Equation  is  an  integrated  expression  arising  from 
Euler’s  analysis  of  fluid  torque  and  mechanical  shaft  work  transferred  between  a  fluid 
and  a  rotor  moving  in  that  fluid.  Since  heat  transfer  is  negligible  in  most  turbomachines, 
the  shaft  work  has  been  equated  to  the  specific  total  enthalpy  transfer  (see  Horlock  [3], 
page  77,  Equation  4.3)  as  follows: 

Ah0  =  A(UVu)  (1) 

where  Vy  is  the  component  of  absolute  velocity,  V,  in  the  direction  U  of  the  rotor  (or 
energy-transferring  device).  An  inconsistency  arises  in  the  consideration  of  a  propeller 
windmilling  on  a  frictionless  shaft  in  a  moving  fluid.  Since  the  shaft  delivers  no  work  a 
viscous  correction  is  necessary  whether  the  right  member  of  (1)  represents  work  or  a 
change  in  the  specific  to'al  enthalpy. 

An  application  of  the  differential  turbomachinery  equation  is  described  for  a  two- 
dimensional,  ideal,  linear  turbine. 


THE  TRANSFORMATION  BETWEEN  ABSOLUTE  AND  MOVING  FRAMES 

In  the  following  discussion  the  subscripts  v  and  w  represent  the  absolute  and  the 
moving  frame  coordinate  and  vector  values.  (See  Figure  1  and  the  Nomenclature  for  defi¬ 
nitions  of  quantities.)  A  frequently  used  relation  connecting  the  absolute  and  moving 
frame  velocities  is 

V=W+U.  (2) 

Following  Spannhske  [4],  spatial  derivatives  in  the  moving  frame  and  time  deriva¬ 
tives  have  the  following  relationship: 


Similar  results  are  obtained  with  cartesian  coordinates. 


Equations  (3)  and  (4)  define,  in  fact,  the  crypto-steady  criterion  [5],  which  if  U  is  con¬ 
stant,  indicates  that  a  frame  exists  in  which  the  flow  regime  may  be  truly  steady  state. 

Since  the  vector  operator  V  is  independent  of  time  and  since  vector  operators  are 
independent  of  frame, 

Vv  =  VH,  =  V  .  (5) 

Now  from  Equation  (3)  for  any  static  function/. 


=  ~U-Vwf=-UVf. 


(6) 


where  q,  represents  all  the  position  coordinates  in  the  absolute  frame  and  is  the  mov- 
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The  Transformation 

ty  —  tyy  =  t 


Fig.  1 .  Configurational  relationships  between  the  absolute  and  moving  coordinate  systems. 
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ing-ftame  coordinate  in  the  direction  uf  the  lifting  body.  Equation  (6)  is  an  extension  of 
the  usual  crypto-steady  relation. 

To  appreciate  the  significance  of  the  final  equation  of  this  section,  considei  an  ob¬ 
server  located  on  a  blade  of  a  rotating  windmill  in  an  infinite  fluid  which  at  infinity 
translates  with  uniform  constant  velocity.  The  observer  cannot  detect  any  time-dependent 
changes  in  the  fluid  at  any  given  point  on  the  blade.  However,  on  moving  to  the  absolute 
frame,  the  observer  notes  changes  in  velocity,  pressure  and  temperature  as  each  blade 
passes  (See  Dean  [6]).  Thus,  the  inequality 


indicates  that  in  the  moving  frame  all  partial  time  derivatives  may  vanish  while  the  abso¬ 
lute  partial  time  derivative  may  be  finite. 

Note,  however,  in  contrast  with  the  inequality  (7)  that  the  substantial  derivative  of 
any  static  quantity  at  a  point  is  independent  of  the  coordinate  system  or  frame  and  is 
therefore  an  invariant  with  respect  to  frame. 

DERIVATION  OF  THE  DIFFERENTIAL  AND  INTEGRAL  FORMS 


UNCOUPLING  THE  SUBSTANTIAL  TOTAL  ENTHALPY  DERIVATIVE  IN 
THE  MOVING  AND  ABSOLUTE  FRAMES 

Relationships  between  absolute  and  relative  flow  fields  are  given  by  the  Coriolis 
form  of  the  Navier-Stokes  Equations  [3,7-10]. 

—  +  W2/2-V  x  (VxD 

dt 


■  +  VW2/2-W  x  (V  x  W)  +  2w  x  W-VU2/ 2 


-Vp  1 _  =, 

—  +  — V -TC'  , 
Q  Q 


where  n'  represents  the  stress  tensor  excluding  the  pressure  tensor,  pd  2  ,  and  U  is 
independent  of  time.  The  relative  acceleration,  dW/dt,  is  defined  in  the  relative  frame, 
i.e.,  (dW/St) w. 

Using  the  absolute-frame  equality  of  (8)  and  the  gradient  form  of  Gibb’s  equation  of 
state,  the  substantial  derivative  of  the  total  enthalpy  is  obtained  in  terms  of  the  partial 
derivative  of  the  pressure,  the  substantial  derivative  of  the  entropic  energy,  T(Ds/Dt) , 


and  the  stress  tensor  n’ .  (See  Wu  et  al.  [7,  8]. 


Dt  Q  \  dt  /  Dt  Q 
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At  this  point  non-ideal  analysis  stops  unless  the  partial  derivative  of  the  pressure 
with  respect  to  time  is  resolved  in  terms  of  the  pressure  gradient  [9].  Employing  Equation 
(6)  with  the  static  pressure  as  the  arbitrary  function  the  result  is 


Dhg 

Dt 


V  „  Ds  V  = 

—  Vp  +  T  —  -4  —  •  (V  ■  7i  )  . 
g  Dt  g 


(10) 


Now  Equation  (8)  may  be  used  to  replace  the  pressure  gradient  in  (10).  From  the  dot 
product  of  U  with  the  moving  frame  equality  in  (8)  the  pressure  gradient  term  becomes 

-  Vp  -  dW  -  —  „  —  _  _ 

-U  •  —  =  (7 - +U-(W ■  V)W  +  2 UaTxW 

g  dt 


-U-VU2/ 2 —  U  ■  (V  n') 


The  gradient  of  the  static  pressure  is  of  course  invariant  in  all  frames,  but  the  mov¬ 
ing  frame  is  convenient  and  preferable  in  (11).  Combining  (10)  with  (1 1),  the  substantial 
derivative  of  the  total  enthalpy  coupled  with  ne  substantial  derivative  of  the  entropic  en¬ 
ergy  is  obtained  in  terms  of  the  moving-frame  flow  field  variables. 


Dh0  Ds  _  aw  _  _ _ _  __  1  _  _  = 

— °--T — =U - +  U  ■  (W  ■  V)W  +  2U  ■  to  x  W  +  -W- (V  n')  . 

Dt  Dr  dt  g 


(12) 


Wu  (reference  [8j,  page  91)  noted  that  in  the  moving  frame  following  a  particle  of 

fluid 


=  q  +  —  (rr'  •  V)  ■  W  , 
Q 


(13) 


where  q  represents  heat  transfer,  and  the  second  term  is  the  specific  dissipation.  If  it  is 
assumed  that  heat  sources  and  sinks,  and  thermal  conduction  are  negligible,  then  from 
(12) and  (13) 

Dh0  _  dW  _  _  „  _  _  _  1  =  _ 

- —  =  U - +{/-(W-  V)W  +  2U  a>  x  W  +  -V  -{rf  W)  .  (14) 

Dt  dt  g 

Equation  (14),  the  essential  development  of  this  paper,  is  a  differential  turbomachin¬ 
ery  equation  expressed  in  terms  of  the  velocity  vector  of  the  moving  frame.  It  is 
universally  applicable  in  any  moving  frame  whether  rotating  or  not. 


THE  SUBSTANTIAL  DERIVATIVE  OF  THE  TOTAL  ROTHALPY 

Using  the  moving  frame  equality  of  Equation  (8)  and  Gibb's  Equation  of  state  for 
the  pressure,  a  viscous  moving-frame  version  of  Crocco’s  Equation  may  be  obtained. 

\  _  _  i  _ 

—  +  VW2/2-W  x  (2m+V  x  W) -  \U2 /2  =  - V/i  +  TVs  +  —  V  ■  rr'  . 

r  e  (15l 
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In  (15)  the  subscript  w  on  the  relative  acceleration  term  implies  that  spatial  vari¬ 
ables  in  the  relative  frame,  r,6w  and  z  are  fixed. 

Wu’s  total  re  ative  rothalpy,  hOWt  may  be  defined  as  follows: 

how  =  h  +  W1/ 2  -  U2/2  =  h0-  UVv  =  h  +  V2/2  -  U\T-j  .  (16) 


Terms  may  be  regrouped  in  (15)  to  obtain  an  expression  explicit  inA/w- 


V/w  =  -f—  I  +W  x  (2aT  +  V  X  R7)  +  TVs  +  —  V  •  zc' .  (17) 

\  dt  j*  e 


Now  the  definition  of  the  substantial  derivative  of  the  total  relative  rothalpy  in  the  mov¬ 
ing  frame  is  as  follows: 


Taking  the  dot  product  of  the  relative  velocity  W  on  (17)  and  adding  the  dot  product  to 
(18)  yields 


+  TW -Vs +  —  ■  {V  -tz') 
Q 


(19) 


Using  Gibb’s  Equation  to  eliminate  the  enthalpy  yields  an  equation  which  is  analo¬ 
gous  to  (9)  in  the  relative  frame. 


W  = 

+  —  (V*')  . 

M'  Q 


(20) 


Since,  as  noted  above,  the  substantial  derivative  of  the  static  quantities  are  invariant 
in  all  frames,  using  (!3)  and  assuming  that  heat  transfer  is  negligible,  (20)  may  be  rewrit¬ 
ten  thus: 


H'  Q 


(21) 


The  quantity  h0„  in  Equation  (21 )  exhibits  an  explicit  time  dependence  in  dp/dt 
which  was  deliberately  ignored  in  crypto-steady  flow  [91  but  recently  reconsidered  [12]. 
Also,  in  (21 )  if  the  viscosity  does  not  vanish,  how  also  depends  upon  the  last  term  [9] 
which  includes  the  viscous  work  and  the  viscous  dissipation. 
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CRYPTO-STEADY  STATE  CONDITIONS 

We  now  examine  the  case  where  the  observer  is  fixed  to  an  infinite  circular  cylinder 
with  bound  circulation  translating  uniformly  in  a  infinite  fluid.  The  observer  fixed  to  the 
moving  cylinder  cannot  sense  any  time  dependence.  Of  course  time-dependent  velocity 
and  pressure  are  immediately  sensed  when  the  observer  moves  to  the  absolute  frame. 

In  another  example  the  observer  is  fixed  to  a  blade  of  an  isolated  rotating  air  screw 
or  a  marine  screw  in  an  infinite  uniform  fluid.  No  time-dependence  can  be  sensed  in  the 
moving  frame  and  a  crypto-steady-state  prevails. 

Under  crypto-steady-state  conditions  the  partial  derivative  of  the  pressure  with  tune 
must  vanish.  This  poses  a  great  simplification  which  justifies  common  practice  [11]  in  the 
design  of  marine  screws  where  time  dependence  is  generally  ignored. 

Before  departing  from  the  crypto-steady-state,  it  is  pertinent  to  ask  what  happens  to 
rotor  energy  transfer  if  there  is  steady-state  in  the  absolute  frame  so  that  dp/ dtv  vanishes. 

Clearly,  the  viscous  terms  in  (9),  (10),  and  (14)  do  not  vanish.  In  (14)  only  the  viscous 
term  of  the  right  member  makes  any  contribution  to  bladeless,  viscously  coupled  discs 
and  concentric  cylinders. 

Now,  if  viscosity  vanishes  and  the  time-dependent  terms  vanish  because  we  impose 
the  crypto-steady  condition,  then  (21 )  reduces  to 

W  ■  Vhow  =  0  .  (22) 

The  non-trivial  solution  is 

V/w  =  V(/1  +  W2/2  -  U2/ 2)  =  V (hQ  -  UVu)  =  0  .  (23 ) 

Equation  (23)  represents  gradient  forms  of  the  classic  Turbomachinery  Equation  in  the 
moving  and  absolute  frame.  Integration  of  the  right  member  of  (23)  over  a  stream  tube  in 
the  absolute  frame  yields  the  classic  Turbomachinery  Equation  (1 ).  The  argument  demon¬ 
strates  that  Equation  ( 1 )  is  strictly  true  only  for  isentropic  and  steady-state  flow. 

CHOOSING  A  FRAME  FOR  TIME  DEPENDENCE 

Since  the  analyses  are  based  upon  two  frames  of  reference,  the  hypothesis  now  pro¬ 
posed  that  time  dependence  may  be  associated  with  disturbances  either  in  the  bulk  flow 
or  with  flow  local  to  the  moving  frame.  Therefore,  formulation  of  a  mathematical  de¬ 
scription  of  the  flow  entering  into  the  calculational  domain  may  be  most  conveniently 
described  in  bulk  frame  coordinates  (the  absolute  frame,  and  subsequently  transferred  to 
the  moving  frame).  On  the  other  hand  flow  disturbances  arising  near  the  moving  body  are 
most  conveniently  described  in  terms  of  the  moving  frame  coordinates.  Since  the  discus¬ 
sion  centers  on  rotors,  it  is  natural  to  consider  subsequent  time  fluctuations  in  the  moving 
rotating  frame. 

A  SINUSOIDAL  VELOCITY  IN  THE  MOVING-FRAME 

Now  in  Equation  (14)  a  small  sinusoidal  moving-frame  oscillation  is  imposed  on  the 
relative  velocity  W  leading  to  a  new  velocity,  W  +  W'  .  The  ratio  of  the  magnitude  of 
time-independent  W”  to  time-independent  W  will  be  arbitrarily  fixed  at  ^  0.05.  The 
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imposed  frequency  of  the  time-dependent  W'  is  at  least  an  order  of  magnitude  greater 
than  the  blade  velocity.  If  the  oscillating  frequency  is  too  low,  it  will  exhibit  a  period  that 
is  long  relative  to  the  blade  period  and  is  the*-* fore  not  considered.  Thus,  prior  to 
Reynolds  time  averaging  (which  should  not  employ  too  large  an  interval  to  avoid  time 
smoothing)  the  velocity  and  other  thermodynamic  quantities  will  fluctuate  at  a  point  (on 
the  blade,  for  example)  during  a  blade  period.  With  Reynolds  averaging,  the  partial  de¬ 
rivative  of  W'  with  respect  to  time  will  vanish. 

However,  the  Reynolds  average  over  non-linear  terms  will  lead  to  non-vanishing 
Reynolds  stress  terms  which  will  not  be  addressed  here.  As  a  result  of  the  Reynolds  time¬ 
averaging  process,  Equation  (14)  may  be  written 
r\L 

-~-  =  U-[(W  V)  W]  +  2U  st  x  {tr}  +  {l/e  V-ar'-IP]  (24) 

where  the  braces  represent  the  Reynolds  time  average  and  W  and  it'  include  the  time- 
dependent  W  .  This  notation  will  not  be  used  again. 

Now,  note  that  arbitrary  samples  of  turbulent  flow  over  short  (relative  to  the  blade) 
periods  may  be  written  as  a  Fourier  summation  of  regular  sinusoidal  oscillations.  For  the 
limited  picture  developed  here.  Equation  (24)  is  again  applicable  by  the  argument  above. 

SOME  CHARACTERISTICS  OF  THE  NEW  DIFFERENTIAL 
TURBOMACHINERY  EQUATION 

Equation  (14),  an  explicit  statement  of  the  substantial  derivative  of  the  total  enthal¬ 
py,  represents  the  power  transfer  at  any  point  in  a  flow  field  between  the  fluid  and  a  bodv 
moving  at  constant  speed  in  that  fluid. 

With  (13)  it  is  now  possible  to  uncouple  the  total -enthalpy,  pressure  relationship  of 
(10)  from  the  entropic  energy. 

^-  =  _^.Vp  +  i-V  (^'  .r)  +  -t7-(V -n')  .  (25) 

Dt  Q  Q  Q 

An  inviscid  Equation  (25)  indicates  that  neither  axial  nor  radial  pressure  gradients 
aie  germane  to  the  calculation  of  specific  total  enthalpy  transfer.  (Mass  flow  rates  are  of 
course  a  function  of  axial  or  radial  pressure  gradients.)  Note,  in  contrast  to  reference  [3], 

p.  7),  for  the  inviscid  case,  only  transverse  pressure  gradients  parallel  to  U  contribute  to 
total  enthalpy  transfer.  The  impulse  stages  of  turbomachines  prove  that  axial  or  radial 
pressure  gradients  need  play  no  role  in  energy  transfer. 

For  bladeless  viscous  coupling  of  discs  or  concentric  cylinders,  the  gradient  of  the 

pressure  parallel  to  U  vanishes.  This  is  equivalent  to  steady-state  pressure  in  the  absolute 
frame  (see  reference  [6]). 

A  proper  test  of  (14)  would  be  whether,  on  integration  over  the  rotor  blade-to-blade 
flow,  it  would  predict  a  total  enthalpy  transfer  compatible  with  that  of  the  classic  Turbo¬ 
machinery  Equation  (1 ).  Therefore,  the  integration  of  (14)  for  crypto-steady  conditions 
will  be  performed  as  a  test  in  the  three-dimensional  domain.  Then  an  application  of  the 
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new  equation  will  be  developed  in  a  two-dimensional  linear  turbine.  Finally,  the  steady- 
state  viscous  problem  of  two  concentric  cylinders  will  be  examined. 

INTEGRATION  OF  THE  TOTAL  ENTHALPY  RATE 

In  the  integration  process  it  will  be  assumed  that  the  flow  may  be  divided  into 
streams  which  pass  between  a  given  pair  of  blades.  In  the  rotating  frame  the  streamtube 
walls  are  fixed  steady-state  walls  associated  with  a  steady-state  mass  flow  rate  m  which 
may  consist  of  radial  and  axial  mass  flow  components. 

DERIVATION  OF  THE  INTEGRAL  FORM  FROM  THE 
DIFFERENTIAL  FORM 

The  differential  form  of  the  turbomachinery  equation  (14)  may  be  integrated  to 

yield 


q  [U  ■  (W-  V)W+  2U- UJ  x  fV]  rdrdQdz 


III 


V  ■  (n'  -W)  rdrdOdz 


The  first  term  of  the  right  member  of  (26)  is  the  tangential  component  of  the  con¬ 
vective  term  obtained  on  dot  multiplication  with  U  ,  i.e.: 
eU(WV)W 


*eU 

\  dr  r  d6  dz  r 


The  first  and  last  terms  of  (27)  will  be  combined  in  an  integral  identified  by  /  ]  4  thus: 


r  c  «/ 

/i,4  =  j  J  j  pray 


■rdrdOdz  . 


The  factors  in  (28)  may  be  rearranged  to  express  the  radial  mass  flow  which  must  be 
constant  in  steady  flow.  Following  arguments  presented  previously  [9] 


=  mr| 


d{UWe)  =  mr[{UWe)r2  -  {UWeU)  , 


where 


Wfl(r)  =  — — —  Wed6dz 

A0Az  J  Jfl 
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The  third  term  of  the  right  member  of  (27)  may  be  written  to  show  the  axial  mass 
flow  rate  m.  explicitly. 


h  = 


gW2rdrd0 


dUWe 

dz 


dz  , 


(31) 


where  the  axial  velocity  Wz ,  V/q  ,  and  U  have  been  averaged  over  r  and  6  .  Substitut¬ 
ing  the  axial  mass  rate  for  the  parenthesis  in  (31), 


=  mz  j 


d(UW0)  =  m2[(UWe)a  -  (t/W0*i]  • 


(32) 


The  second  term  of  the  right  member  of  (27)  provides  an  integral  which  contains  the 
tangential  kinetic  energy. 


^-dQ  \drdz 
30 


■ll‘r(l 


=  I  [w»(02) "  wi{6^]drdz  =  0 


(33) 


Since  the  tangential  velocities  at  the  blade  walls  are  equal  to  the  blade  velocity,  the 
integral  vanishes. 

Now  identifying  the  second  term  of  the  right  member  of  Equation  (27)  as  I5,  we 
may  write 


=  11/ 


2garrWrrdrd0dz  . 


Finally,  following  the  arguments  above, 

h  »  Wr[(f/2)r2-(^2)rl]  • 

Summing  the  components  of  integration  1 14  through  1$, 


(34) 


(35) 


\\  \Q  1**  =  mr  Ar(U^e)  +  AAU2) 


+  m2A2(UWe)  +  (VIS) 


(36) 


where  Ar  and  Az  represent  the  change  along  r  and  z  respectively,  and  (VIS)  is  the  inte¬ 
grated  viscous  term. 

Now  adding  mzAz()2,  which  is  zero,  to  Equation  (37), 
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(37) 


f  f  |e^-rfr  =  m^[(f/  +  ^)C/] 

+  mzA2[(f/  +  ]&g)U ]  +  (V7S)  . 

In  Equation  (37)  the  terms  Wg  are  averaged  over  6  and  z  in  the  first  term  and  over 

r  and  6  in  the  second  term.  If  the  total  steady-state  mass  rate  m  between  a  pair  of  blades 
is 

m  -  mr  +  mz  ,  (38) 

then, 

A ho  =  frbrdrie)  +fz&z(lrie)  +  ( VIS)/m  .  (39) 

Equation  (39)  represents  classic  Turbomachinery  Equation  with  mixed  flows  and 

the  terms  Vg  and  U  are  averaged  over  the  blade  space  where  necessary.  The  coeffi¬ 
cients  £  and  fz  represent  the  radial  and  axial  fractions  of  the  mass  flow. 

The  viscous  term  is  a  novel  feature  of  the  derivation  which  explains  losses  in  rotors 
during  windmilling.  With  the  exception  of  the  viscous  energy  term  (VIS),  the  integral  ex¬ 
pression  equation  (39)  exhibits  a  formal  similarity  and  compatibility  with  the  classic 
Turbomachinery  Equation.  The  differences  arise  because  the  integration  in  (39)  has  been 
performed  over  the  entire  blade-to-blade  volume  rather  than  a  stream  tube.  It  is  easy  to 
combine  the  axial  and  radial  flow  terms  by  letting  A  vary  in  both  r  and  z  .  The  deriva¬ 
tion  lends  credence  to  the  hypothesis  that  Equation  (14)  is  indeed  a  differential 
turbomachinery  equation. 

A  two-dimensional  application  and  test  of  the  differential  form  (14)  on  an  ideal  lin¬ 
ear  device  where  the  solution  is  known  precisely  will  now  be  examined. 

THE  SUBSTANTIAL  TOTAL  ENTHALPY  RATE  IN  A 
TWO-DIMENSIONAL  DEVICE 

An  infinite  circular  cylinder  with  bound  circulation,  as  shown  in  Figure  2,  is  an  ele¬ 
mental  linear  turbine.  It  may  be  considered  as  an  infinite  sail  on  a  sailboat  or  an  infinite 
wing  on  a  sailplane.  The  device  extracts  energy  from  the  ideal  inviscid  working  fluid. 
Work  is  performed  on  the  sailplane  (fixed  to  a  vertical  rail)  by  raising  its  height  at  uni¬ 
form  speed  U  against  gravity.  Work  on  the  sailboat  is  performed  by  moving  the  boat  at 
uniform  speed  U  which  elevates  a  weight  attached  at  minus  infinity  by  an  infinite  tether. 
In  the  moving  frame  the  apparent  velocity  of  the  ideal  working  fluid  at  infinite  distance  is 
W0  .  The  relationship  between  the  absolute  and  moving  coordinate  system  and  the  velo¬ 
cities  is  given  by  the  transformation  of  Figure  1 . 

Since  the  flow  field  is  ideal,  the  flow  domain  may  be  described  by  a  potential  func¬ 
tion  or  its  conjugate  stream  function.  The  lift  is  therefore  the  ideal  lifting  force,  L,  of  the 
Kutta-Joukowski  Equation  given  by 
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L  -  qW0T  , 


(40) 


where  T  is  the  scalar  circulation.  The  units  are  force  per  unit  length  of  cylinder.  In  the 
absolute  and  moving  frame  the  lift  component  Ly  directed  parallel  to  the  y  axis  of  Figure 
2  is  given  by 

Ly  =  QWoxr  =  qVoxT  ,  (41) 

where  the  subscript  x  represents  the  x  component.  Recalling  that  U  is  the  velocity  of  mo¬ 
tion  of  the  device  (sail  or  wing  or  rotating  cylinder)  as  perceived  in  the  absolute  frame, 
the  power  is  the  product  of  U  and  Ly. 

Power/ unit  length  =  qUWoxT  .  (42) 

Since  we  assume  that  there  is  no  heat  rate, 

DH 

--JL  =  ^qUW„xT  ,  (43) 

where  H0  is  the  total  enthalpy  of  the  system  per  unit  length.  Equation  (43)  is  the  antici¬ 
pated  relationship  which  should  ultimately  be  developed  from  the  differential  form  (14). 

THE  STREAM  FUNCTION,  VELOCITY  AND  RELATIVE  ENTHALPY 
IN  THE  FRAME  OF  THE  BLADE 

Since  ideal  flow  has  been  assumed  in  the  moving  frame  of  the  blade,  the  stream 
function,  ^ ,  is  the  usual  function  modified  for  motion  along  the  y  axis. 

V'  =  -  Woy(  1  -  a2/ r^)r  cos  6  +  Wox(  1  -  a2 /rV  sin  0 


+  (T/2ji)ln(r/a)  . 


(44) 


The  constant  a  is  the  radius  of  the  cylinder.  The  cartesian  velocity  components  are 
obtained  by  the  usual  transformation  [9]  as  follows: 


Wx  = 


Wy  = 


w  ,  a2W0X(y1-x2)  2a2Woyxy 
(x2  +  y2)2  (x1  +  y2)2 

+  JL  _ 1 _ 

2nr  (x^  +  y2)  ’ 

la2Woxxy  aty^-y2) 
^  (x 2  +  y2)2  ( x 2  +  v2)2 

_T_  x 

2 it  {x 2  +  y2) 


(45) 


(46) 


Now  the  relative  vorticity  must  vanish  because  potential  flow  cannot  have  vorticity. 
A  check  of  the  vorticity  in  the  relative  frame  shows  that  indeed  it  vanishes.  Also,  the 
time-dependent  term  vanishes. 
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THE  SUBSTANTIAL  TOTAL  ENTHALPY  DERIVATIVE 

In  the  linear  two-dimensional  system,  the  differential  form  of  the  turbomachinery 
Equation  (14)  is  simplified  because  the  rotation  vanishes. 

—  =U(WV)W.  (47) 

Dt 


Since  the  vorticity  vanishes 


W  VW  =  VW2/2  , 


(48) 


and 


—  =  U-  VW2/2  = 
Dt 


U  dW2 

2  dyw 


(49) 


The  integrated  substantial  total  enthalpy  rate  per  unit  length  [where  the  subscript  on  y  in 
(49)  has  been  dropped]  is 


DH0 

Dt 


q£  d(W2+W2) 
2  dy 


dydx  . 


(50) 


Integration  of  (50)  will  be  performed  over  all  space  per  unit  length  z  of  the  blade.  The 
choice  of  time  is  immaterial  since  the  fluid  dynamics  are  steady  state  in  the  moving  frame 
and  the  thermodynamic  rates  over  all  space  are  invariant  with  time.  It  is  understood  that 
the  integration  applies  only  to  the  fluid  domain  and  that  boundaries  at  solid  walls  are  ob¬ 
served. 


Now  without  going  into  details  [9]  the  integral  of  the  total  derivative  in  (50)  is 


DH0 

Dt 


=  -QU 


4UWox  x(a2  -  x2)1/2(2x2  -  a2) 


2WoxTx{a2  -  x2)1/2  WoxUx{a2  -x2)1'2 

+  - A -  + - 7 - 


Wox rV-x2)1/2  ,  4W0XUxia2-x2)1/\a2-lx2) 
+  +  .  -  —  - 


WoxT(a2-x2)l/2{a2-2x2)\j 

+ - ^ - ]dx 

Jia  I 

Note  that  only  odd  terms  in  y  make  any  contribution  to  (51).  Since  the  first  and  fifth 
terms  cancel,  only  four  terms  remain.  The  integration  with  respect  to  x  is  performed 
through  a  transformation  employing 
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x  =  a  cos  0  , 


(52) 


with  integration  limits  given  by 

6  =  n  when  x  =  -  a  , 
0  =  0  when  x  =  a  . 

Making  the  substitutions 


(53) 


1° 


8 WoxUa  cos  6  sin2  6d6  - 


Wo£  f° 


sin2  BdO 


WoxT 

JT 


/: 


sin2  0(1  -  2  cos2  6)dQ  1  . 


(54) 


In  (54)  the  second  integral  makes  no  contribution  because  it  is  antisymmetric.  The  first 
integral  cancels  the  second  term  in  the  last  integral  to  yield  from  the  surviving  terms 

DH 

—r 1  =  -  qW0XYU  =  -  gVoxTU  .  (55) 

Dt 

Equation  (55)  is  identical  with  (43)  and  this  result  illustrates  a  useful  application  of 
the  differential  form  and  constitutes  confirmation  of  the  validity  of  the  differential  turbo¬ 
machinery  Equation  (14).  For  the  linear  case,  the  energy  transfer  rate  of  the  rotor  is 
proportional  to  the  component  of  the  kinetic  energy  gradient  parallel  to  the  moving  rotor. 
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new  equation  will  be  developed  in  a  two-dimensional  linear  turbine.  Finally,  the  ateady- 
state  viacous  problem  of  two  concentric  cylinders  will  be  examined. 

INTEGRATION  OF  THE  TOTAL  ENTHALPY  RATE 

In  the  integration  process  it  will  be  assumed  that  the  flow  may  be  divided  into 
streams  which  pass  between  a  given  pah  of  blades.  In  the  rotating  frame  the  streamtube 
walls  are  fixed  steady-state  walls  associated  with  a  steady -state  mass  flow  rate  m  which 
may  consist  of  radial  and  axial  mass  flow  components. 

DERIVATION  OF  THE  INTEGRAL  FORM  FROM  THE 
DIFFERENTIAL  FORM 

The  differential  form  of  the  turbomachinery  equation  (14)  may  be  integrated  to 

yield 


V)tP+  2U-  fiT  x  W\  rdrdOdz 


t 

\ 


III 


V  ■(jf'  -W)  rdrdOdz 


The  first  term  of  the  right  member  of  (26)  is  the  tangential  component  of  the  con¬ 
vective  term  obtained  on  dot  multiplication  with  V ,  i-e.: 

qU(WV)W 

\  lr  r  d$  »z  r  f  Q 

The  first  and  last  terms  of  (27)  will  be  combined  in  an  integral  identified  by  / 1,4  thus: 

'1.4  -  /  /  Je  m  ^  Jp-rdrjedi  (2 

mf(fl*r'm*)£l?L*- 


where  we  have  used  the  mean  value  theorem  to  take  W(r)  outside  the  double  integral. 
Thus, 

^  I  C  (29> 

The  fartor  in  parentheses  in  the  right  member  of  (28)  is  the  radial  mass  flow  at  any  point 
r. 
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"Hr)  ~  fl  QWfTdOdz  =  mfr(r)  , 

where  Mr)  is  the  radial  fraction  of  tne  mass  flow  rate  m.  Then 

'  hA  ««  J Mr )~p-  dr  =  m  j  fJ(U&e)  . 

Now  in  the  steady  state  the  nuu  w*Iue  of  Mr) ,  fr  ,  is  a  constant  given  by 

fexk  _ 

f,WWe)d(UWs)/A(inir:  \ 

'  ]  menace 

and. 


/ l,a  *  mfrMUWg)  =  -  (CW,)r,] 


(30) 


(31) 


(32) 


(33) 


Using  the  same  arguments  as  used  for  the  radial  mass  flow,  the  third  term  of  the 
right  member  of  (27)  may  be  written  to  show  the  axial  mass  flow  rate  m,(r)  explicitly. 

(See  argument  above  for  mAr) .)  With  die  velocities  Wz ,  Wq  ,  and  (7  averaged  over  r 
and  6  ,  one  obtains 

/3-m*[(^«)%-(C^)r,J  .  (34) 

The  second  term  of  die  right  member  of  (27)  provides  an  integral,  l , ,  which  coo- 
tains  the  tangential  kinetic  energy. 


- / .  (35) 

Since  die  tangential  velocities  at  the  blade  walls  are  equal  to  the  blade  velocity,  the 
integral  vanishes.  f 

Now  identifying  the  second  term  of  the  right  member  of  Equation  (27)  as  Is,  we 
may  write 

/3  -  J  J  J  iQwb-Wrrdrdedz  *  mr[(U\  -  (U\)  .  (36) 

Summing  the  components  of  integration  (33)  through  (36) 

Jjj  + 

*191^0^9)  + (VIS)  .  (37) 
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where  Ar  and  A2  represent  the  change  along  r  and  z  respectively,  and  (VIS)  is  the  inte¬ 
grated  viscous  term. 

Now  adding  mzA:U2,  which  is  zero,  to  Equation  (37), 

IJ 


( U+We)U 

In  Equation  (3&)  the  terms  Wg  are  averaged  over  6  and  z  in  the  first  term  and  over 
rand  6  in  the  second'  term.  Finally, 

the  +/A*(UVe)  +  (VIS)/m  =  MWe)  +  (VIS)/m  .  (39) 

Equation  (39)  repre^uts  the  classic  Turbomachinery  Equation  with  mixed  flows 

and  the  terms  Vg  and  U  are  avenged  over  the  blade  space  where  necessary. 

The  viscous  term  is  a  novel  feature  of  the  derivation  which  explains  losses  in  rotors 
during  windmilling.  With  the  exception  of  the  viscous  energy  term  (VIS),  the  integral  ex¬ 
pression  equation  (39)  exhibits  a  formal  similarity  and  compatibility  with  the  classic 
Itirbomachinety  Equation.  The  differences  arise  because  the  integration  in  (39)  has  been 
performed  over  the  entire  blade- to-biade  volume  rather  than  a  stream  tube.  It  is  easy  to 
combine  the  axial  and  radial  flow  terms  by  letting  A  vary  in  both  r  and  z .  The  deriva¬ 
tion  lends  credence  to  die  hypothesis  that  Equation  (14)  is  indeed  a  differential 
turbomachinery  equation. 

A  two-dimensional  application  and  test  of  the  differential  form  (14)  cm  an  ideal  lin¬ 
ear  device  where  the  solution  is  known  precisely  will  now  be  examined. 

THE  SUBSTANTIAL  TOTAL  ENTHALPY  RATE  IN  A 
TWO-DIMENSIONAL  DEVICE 

An  infinite  circular  cylinder  with  bound  circulation,  as  shown  in  Figure  2,  is  an  ele¬ 
mental  linear  turbine.  It  may  be  considered  as  an  infinite  sail  on  a  sailboat  or  an  infinite 
wing  on  a  sailplane.  The  device  extracts  energy  from  the  ideal  inviscid  working  fluid,  y  • 
Work  is  performed  on  the  sailplane  (fixed  to  a  vertical  rail)  by  raising  its  height  at  uni¬ 
form  speed  U  against  gravity.  Work  on  the  sailboat  is  performed  by  moving  the  boat  at 
uniform  speed  U  which  elevates  a  weight  attached  at  minus  infinity  by  an  infinite  tether. 

In  the  moving  frame  the  apparent  velocity  of  the  ideal  working  fluid  at  infinite  distance  is 
W0 .  The  relationship  between  the  absolute  and  moving  coordinate  system  and  the  velo¬ 
cities  is  given  by  the  transformation  of  Figure  1. 

Since  the  flow  field  is  ideal,  the  flow  domain  may  be  described  by  a  potential  func¬ 
tion  or  its  conjugate  stream  function.  The  lift  is  therefore  the  ideal  lifting  force,  L,  of  the 
Kutta-Joukowski  Equation  given  by 


+  m,A: 


1 


+  (VIS)  =  mAho 


(38) 
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